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ABSTRACT 


Almost  a decade  ago  Sherbrooke  formulated  the  well  known  METRIC 
model  for  determining  optimal  stock  levels  for  recoverable  items  for  two 
echelon  inventory  systems  [3].  Subsequently  Fox  and  Landi  [2]  proposed  a 
Lagrangian  approach  for  obtaining  item  stock  levels  for  each  location. 

In  this  paper  we  develop  a method  for  estimating  the  value  of  the  optimal 
Lagrangian  multipler  used  in  the  Fox-Landi  algorithm,  present  alternative 
ways  for  determining  system  stock  levels,  and  compare  these  proposed 
approaches  with  the  Fox-Landi  algorithm  and  other  solution  techniques.  The 
conclusion  of  this  study  is  that  the  proposed  approximation  methods 
significantly  reduce  computation  time  for  determining  system  stock  levels 
without  degrading  the  quality  of  the  solution. 


r 


z' 


1.  INTRODUCTION 


Almost  a decade  ago  Sherbrooke  formulated  the  well-known  METRIC  model 
for  determining  optimal  stock  levels  for  recoverable  items — items  subject 
to  repair  when  they  fail — in  a two-echelon  setting  [3].  In  particular, 
he  studied  the  Air  Force's  two-echelon  supply  system.  This  system  consists  of 
a set  of  bases  and  their  supporting  depots.  Primary  demands  occur  at  the 
bases  while  depots  are  central  repair  and  inventory  stocking  points  which 
resupply  bases  when  necessary.  When  a failure  occurs  at  a base,  a demand  is 
placed  on  the  base  supply  organization  for  a corresponding  replacement  part. 
Depending  on  the  nature  of  the  failure , the  failed  part  is  then  either  repaired 
at  that  bane  or  is  sent  to  a depot  for  repair.  Resupply  of  the  base  supply 
organization  comes  from  the  base  maintenance  organization  if  repair  is 
accomplished  at  the  base;  othenrise,  resupply  comes  from  a depot.  In  either 
case,  the  organization  resupplying  the  base  supply  activity  does  so  by 
exchanging  on  a one-for-one  basis  a serviceable  part  for  the  failed  part.  Thus 

the  inventory  policy  for  placing  orders  on  the  base's  maintenance  organization 
or  a depot  is  an  (s-l,s)  policy. 

Sherbrooke  presented  a model  (METRIC)  for  determining  both  depot  and 
base  stock  levels  for  all  items  for  this  system.  In  particular,  the  problem 
he  formulated  was  to  minimize  the  average  total  number  of  base  backorders 
existing  at  an  arbitrary  point  in  time  subject  to  a constraint  on  system 
investment;  that  is. 


m n 

minimize  111  (x-s . . )p(x| X . .T . . (s  ) ) 

-L  -L  L i]  ' i]  i]  lO 


j = 1 i=l  x>s. . 

il 


(1) 


m n 


l l c.s  . < C, 
j=0  i=l  1 13 
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subject  to 
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where  n 


represents  the  number  of  items, 


m represents  the  number  of  bases, 

s^_.  represents  the  stock  level  at  base  j for  item  i, 

represents  the  depot  stock  level  for  item  i, 
represents  the  expected  daily  demand  rate  for  item  i at 
base  j , 

c^  represents  the  unit  cost  for  item  i, 

C represents  the  budget  constraint, 

T„(si0)  represents  the  average  resupply  time  for  base  j for  item  i 
given  the  depot  stock  level  for  item  i is  s^Q, 
and  p(x|y)  represents  the  probability  that  x units  are  in  the  resupply 
system  given  that  the  expected  number  of  units  in  the  resupply 
system  is  y. 


Furthermore,  Sherbrooke  shows  that  T. .(s._)  can  be  expressed  as 

13  1O 


T . . (s . ) = r. .A. . + (1-r. . )(B. . + 6(s  )-D.), 

ID  i0  13  13  13  13  i0  1 


where  A„  is  the  average  base  repair  time  for  item  i at  base  j , 

r^.  if  the  proportion  of  demands  requiring  base  repair  for  item  i 
at  base  j , 

is  the  average  depot  to  base  order-and-ship  time  at  base  j for 
item  i , 

is  the  average  depot  repair  cycle  time  for  item  i, 

6(s.q)*D.  = j-  l (x-s^pCxIl.D^,  the  expected  delay  per  depot 
1 1 1 X>SiO 


demand  for  item  i, 
m 

and  X.  = £ (l-r..)X..,  the  expected  daily  depot  demand  rate  for  item  i. 

1 j=l  13  13 


In  the  remainder  of  the  paper  i will  refer  to  an  item  and  j will  refer  to 
a base  (j=0  represents  the  depot);  thus  i and  j will  always  be  elements 
of  the  sets  {l,...,n}  and  {0,...,m},  respectively.  Additionally,  an 
integer  k appearing  in  the  text  to  the  right  of  the  statement  of  a problem 
or  equation  will  designate  for  future  ref erence  that  problem  or  equation.  For  a 
complete  description  of  this  problem's  background  and  formulation  see  reference 
3. 

Subsequently  Fox  and  Land!  suggested  a Lagrangian  approach  for  solving 
Problem  1 [2].  One  of  the  major  obstacles  to  the  successful  implementation  of 
METRIC  using  the  Fox-Landi  algorithm  is  the  requirement  of  estimating  an 
appropriate  value  for  the  Lagrangian  multiplier.  An  important  and  related 
problem  is  the  lengthy  computer  run  time  required  to  obtain  an  optimal  solution 
to  Problem  1 when  using  this  algorithm. 

The  purposes  of  this  paper  are  to  present  an  approach  for  obtaining  an 
estimate  of  the  optimal  Lagrange  multiplier  value  required  in  the  Fox-Landi 
algorithm,  to  present  two  new  methods  for  determining  stock  levels,  and  to 
compare  these  methods  with  the  Fox-Landi  method  and  other  techniques.  The 
proposed  approach  eliminates  the  particularly  time  consuming  portion  of  the 
Fox-Landi  algorithm  devoted  to  searching  for  the  best  Lagrange  multiplier  value. 
Ihe  conclusion  of  the  study  is  that  the  proposed  approximation  methods  signifi- 
cantly reduce  computation  time  for  determining  stock  levels  without  degrading 
the  quality  of  the  solution. 


2.  THE  APPROXIMATION  PROBLEM 


We  begin  this  section  by  constructing  a problem  that  is  a continuous 
approximation  to  Problem  1.  We  next  state  and  prove  two  theorems  that  are 
the  basis  for  an  algorithm  developed  in  the  next  section  which  can  be  used 
to  solve  this  approximating  problem. 


Recall  that  the  total  average  base  backorders  existing  at  any  point  in 


time  for  item  i can  be  expressed  as 


m 

1JS  <*-»!,  >P<xl»ljTiJ(sI0)). 

ij 

Two  useful  probability  distributions  for  describing  the  demand  process  are  the 

Poisson  and  negative  binomial  distributions.  As  shown  in  reference  1,  this 

implies  that  if  demand  has  a Poisson  or  negative  binomial  distribution,  then 

for  a given  X..  T..(s.„),  the  form  of  p(x| A . .T. . (s. - ) ) , the  probability 
l]  l]  lO  r ' ij  i]  iO  r J 

distribution  representing  the  number  of  units  in  resupply  of  item  i at  base 
j at  any  point  in  time,  is  a Poisson  or  negative  binomial  distribution, 
respectively. 

Experimental  data  gathered  during  the  conduct  of  this  study  indicate 
that  when  p(x| ^sio^  ^ e-’t*ier  a Poisson  or  negative  binomial 
distribution,  the  above  total  expected  backorder  expression  can  be  closely 
approximated  by  an  exponential  function.  That  an  exponential  function  accurately 
approximates  this  expression  should  not  be  entirely  unexpected.  First,  for 
budgets  of  practical  interest  the  item  stock  levels , s^ , are  normally  much 
larger  than  the  average  demand  during  the  resupply  time.  For  example,  the 
probability  of  running  out  of  stock  during  the  resupply  time  is  often  much  less 
then  .15  in  real  Air  Force  applications.  Thus,  the  only  probabilities  entering 
the  backorder  calculation  are  the  tail  probabilities.  In  the  tails,  the 
Poisson  and  negative  binomial  distributions  behave  almost  like  the  geometric 
distribution,  that  is,  each  succeeding  probability  is  roughly  a constant 
proportion  of  its  predecessor.  Consequently,  when  s^.  is  large  relative 
to  A (s^0>,  the  expected  number  of  backorders  existing  at  any  point 
in  time  at  location  j for  item  i is  approximately  a geometric  function  of  s^.. . 
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Recall  that  the  total  average  base  backorders  existing  at  any  point  in 
time  for  item  i can  be  expressed  as 


V J (x-s. . )p(x| X . .T. . (s.  )) 
A L n r 1 lj  lj  lO 


j=l  x>s.. 
i] 


Two  useful  probability  distributions  for  describing  the  demand  process  are  the 
Poisson  and  negative  binomial  distributions.  As  shown  in  reference  1,  this 
implies  that  if  demand  has  a Poisson  or  negative  binomial  distribution,  then 
for  a given  X . . T..(s.-),  the  form  of  p(x|  X . .T. . (s . ) ) , the  probability 
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distribution  representing  the  number  of  units  in  resupply  of  item  i at  base 
j at  any  point  in  time,  is  a Poisson  or  negative  binomial  distribution, 
respectively. 

Experimental  data  gathered  during  the  conduct  of  this  study  indicate 
that  when  p(x|x^T„  (s^Q))  is  either  a Poisson  or  negative  binomial 
distribution,  the  above  total  expected  backorder  expression  can  be  closely 
approximated  by  an  exponential  function.  That  an  exponential  function  accurately 
approximates  this  expression  should  not  be  entirely  unexpected.  First,  for 
budgets  of  practical  interest  the  item  stock  levels,  , are  normally  much 

larger  than  the  average  demand  during  the  resupply  time.  For  example,  the 
probability  of  running  out  of  stock  during  the  resupply  time  is  often  much  less 
then  .15  in  real  Air  Force  applications.  Thus,  the  only  probabilities  entering 
the  backorder  calculation  are  the  tail  probabilities.  In  the  tails,  the 
Poisson  and  negative  binomial  distributions  behave  almost  like  the  geometric 
distribution,  that  is,  each  succeeding  probability  is  roughly  a constant 
proportion  of  its  predecessor.  Consequently,  when  s^  is  large  relative 
to  A.jT.j(s.Q),  the  expected  number  of  backorders  existing  at  any  point 
in  time  at  location  j for  item  i is  approximately  a geometric  function  of  s.^ . 


Therefore  an  exponential  function  is  a useful  continuous  fipproximation  to 
this  relationship  between  expected  backorders  at  a location  and  the  item's 
stock  level  at  that  location. 


Furthermore,  total  expected  base  backorders  exhibit  this  same  behavior. 

If  demand  has  either  a Poisson  cr  negative  binomial  distribution  (or,  for 
that  matter,  any  other  compound  Poisson  distribution),  then  the  total  number 
of  units  of  an  item  in  resupply  across  all  bases  also  has  a Poisson  or  negative 
binomial  distribution,  respectively,  given  independence  of  demand  among  bases 
and  assuming  the  "order  size"  distribution  is  the  same  at  all  bases.  Since 
in  most  practical  situations  total  system  stock  substantially  exceeds  the 
total  expected  number  of  units  in  resupply,  the  tail  of  the  distribution 
describing  the  total  number  of  units  in  resupply  is  the  only  portion  of  the 
distribution  of  importance.  As  an  approximation,  this  distribution  can  be 
used  to  determine  the  nature  of  the  relationship  between  total  expected  base 
backorders  and  total  system  stock.  For  the  reasons  discussed  previously,  an 
exponential  function  should  adequately  represent  this  relationship  as  well. 

Thus  we  will  approximate  total  system  backorders  for  item  i,  that  is. 


l l 


i=i  x>s . . 

11 


with  an  exponential  function  of  the  form 


-b.N. 

B . ( N . ) = a . e 1 1 . 
ii  i 


In  this  approximation  IL  represents  total  system  stock.  The  parameters 
a^  > 0 and  b^  > 0 are  estimated  using  regression  analysis.  The  data  used 
in  the  regression  analysis  are  the  backorder  data  obtained  fron  the  solution  to 
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the  problem 


minimize  V V (x-s..)p(x|X..T..(s.  )) 
j=l  x>g  13  13  13  1 

ij 


subject  to 


m 


Y s . . = N . , and 

jio  ^ 1 


s^_.  — Ojl)»  • • )N,, 


for  several  appropriate  values  of  N . 

We  now  formulate  a continuous  approximation  to  Problem  1 in  which  the 
exponential  representation  of  total  system  backorders  for  an  item  is  used. 

In  this  approximation  problem  the  decision  variables  are  the  total  system 
stock,  h\,  rather  than  the  stock  levels  for  each  location,  s „ . As  we  shall 
see,  the  main  reason  for  studying  this  approximation  problem  is  that  it  is 
a vehicle  for  obtaining  an  estimate  of  the  optimal  Lagrangian  multiplier 
value  used  in  the  Fox-Landi  algorithm.  The  approximation  problem  is 
formulated  as 


n 

minimize  / B.(N.) 

i=l  1 1 


subject  to 


I 

i=l 


c.N.  < C, 
11 


N.  > 0. 

1 - 


(2) 
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f 

t 

If 

fc 

f 


I... 

i I 


k 


’t  > 4 


i'J 
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Pm 


F ■ 

T.i 

i A 


Note  that  is  a continuous  variable  in  this  approximation.  The  optimality 


conditions  (Kuhn-Tucker  conditions)  for  this  problem  are  as  follows: 


Find  9 ^ ^ 0 such  that 


dB. 

a) 

1 

dN. 

l 

+ ei' 

n 

b) 

l 

c.N. 

i=l 

n. 

> o. 

c)  9 ( l c.N.  -0=0, 
i=l  1 1 


dB. 


and 


d)  N.(~  + 0.c.)  = 0. 
x dH.  1 l 

l 


A relaxed  version  of  Problem  2 in  which  the  non-negativity  constraint 
on  the  item  stock  level  is  removed  is 


linimize  7 B.(N.) 

i=l  1 1 


subject  to 


(3) 


y c.N.  < C. 
L 11“ 
i=l 


The  optimality  conditions  for  this  problem  are: 


Find  £*2  * 0 such  that 


dB. 

a) 

l 

dH. 

l 

n 

b) 

I< 

i=l 

c) 

92  ( 

+ = °> 


= 0, 


i=l 


! ->  - 
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d)  IK(-r-f  + 0 c.)  = 0. 

1 V 2 1 


We  now  explore  the  relationship  between  Problems  2 and  3 in  detail. 

Suppose  we  obtained  a solution  to  Problem  3 (we'll  show  how  to  find 

its  solution  in  the  next  section).  Let  IK  represent  the  optimal  solution 

2 

to  Problem  2,  and  IK  represent  the  optimal  solution  to  Problem  3.  If 
2 12 

IK  10  for  all  i,  then  IK  = IK  and  the  objective  function  values  are 
i li 


equal. 


Suppose,  however,  that  N.  < 0 for  at  least  one  value  of  i.  Let 


N.  = max(0,N. ) 
l i 


c = 1 c.N. 

i=l  1 1 


— 2 — 2 — 

Since  IK  > N.  for  all  i and  IK  > N.  for  at  least  one  i,  C > C. 

i l li 

Suppose  Problem  2 is  modified  slightly  so  that  the  right-hand  side 
value  C is  replaced  by  C.  This  modified  problem  is 


minimize  T B.(N.) 

i=l  1 1 


subject  to 


l c;Ni  .4  c> 
i=l  1 1 


IK  > 0. 

l • 


The  optimality  conditions  for  this  problem  are  the  same  as  those  given 
for  Problem  2 after  substituting  C for  C.  Also,  let  0 represent  the 


optimal  value  of  the  Lagrangian  multiplier  for  Problem  4. 


In  solving  Problen  3,  we  will  cl  tain  a value  for  62*  We  now  show  that 

— o 

6 = ©2^  and  that  = max(0,N^)  is  the  optimal  solution  to  Problem  4 by 
demonstrating  that  these  values  satisfy  the  Y ihn-Tucker  conditions 
corresponding  to  Problem  4. 


By  construction, 


l c.N  = C,  N.  >0,  and  e(  7 c.N.  - 

i=l  i i i - ili  1 x 


C)  = 0. 


— — — ? 2 
If  6 = 6 , 0 > 0 since  0„  > 0.  Suppose  N.  = N.;  that  is,  II.  >0.  The 

2 2 — ^ li  l — 


dB. 

•—  dlii  2 

N.-N.  N.=N7 

ii  li 


1 M.=H? 

1 l l 


dB.  j 

+ °2Ci  = diTl 


+ 0c. 

i 


N.=N. 
1 i l 


By  assumption  there  exists  at  least  one  value  of  i for  which  lit  > 


that  is,  N.  = 0 while  lit  < 0.  Since 
i i 


dW.  dll.  0 

1 II. =0  1 1 II.  =11? 

1 i 'ii 


due  to  the  exponential  form  of  B^OO,  and 


dB.  d3i 

-r— + 0„c.  =0,  we  know  that  -jrr- 
dh.  _ 2 l dil . 

1 N.=N7  N . = 0 

ii  i 


+ 0c.  > 0. 
i 


Consequently,  the  optimal  solution  to  Problem  4 is  IL  = Ih  = max{0,51.} 
Furthermore,  the  optimality  conditions  are  satisfied  when  0 is  equal  to  0 


r 
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Theorem  1.  0^  >_  0^. 

Proof:  The  optimal  objective  function  value  for  Problem  2 is  a convex, 

differentiable,  strictly  decreasing  function  of  the  available  budget,  C. 

Since  the  slope  of  this  function  is  equal  to  the  negative  of  the  Lagrangian 
multiplier  value,  0^  > 0 since  C < C.  But  q2  = 0,  so  0^  ©2 - | 

Corollary.  0^  > 0^  when  C > C. 

Next  we  compare  IId  with  II..  If  C = C,  then  Nd  = N.  for  all  i. 

^ i l ii 

New  let  us  suppose  C > C so  that  0^  > 0^  = B.  bet  us  examine  the  two  cases 
N^  > 0 and  IK  ~ 0 separately. 

First  assume  IK  > 0.  Then 


v 


!*  • 
' J 


* 

H 


i 


if 

I 


Ll 


dBil 

dir!  + eci = °- 

1 |n.=ii. 


dB. 


Furthermore,  if  II.  > 0,  then  377— 

1 U 14  • 


N.-N^ 
1 1 


+ 0,  c . = 0 . 
1 1 


hB 


Since  0,c.  > 0c.  = - ... 

1 1 1 dll . 


dB. 

dB. 

1 

1 

’ dM . 

> dliT 

II  .=11.  1 

II.  =N.  1 

1 1 


1 1 


IK=Nd 


, and  IId  < N. 
’ 1 1 


1 1 


If  N*  = 0,  then  II.  > Nd. 
1 ’ 11 


"ext  assume  IK  = 0.  Since 


^i 


dll. 

1 


jIK=0 


dB. 

+ 0C.  T-^' 

1 1 di.. 

1 


+ Cc.  > 0,  it  follows  that  ll.  = 0 
1 1 


N.=0 

1 


J 


by  CGT.niOT.Qr.taFy  slackness*  Tvi.'3  we  have  proven  the  following  theere 
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Theorem  2.  fL  ? add it . 1 orally , > whenever  N^  > 0. 

In  this  section  we  established  several  important  relationships  among 
Problems  2,  3,  and  4.  In  the  next  section  we  develop  a simple  algorithm  for 
solving  problem  2 based  on  these  relationships.  We  will  begin  the  next 
section  by  showing  how  to  find  the  solution  to  Problem  3.  As  we  have  just 
demonstrated,  once  we  have  the  solution  to  Problem  3 we  also  have  the  solution 
to  Problem  4.  From  Theorem  2,  we  then  have  an  upper  bound  on  the  value  of 
N^.  In  particular,  if  IP  = 0,  then  = 0.  Combining  this  observation  with 
the  implications  of  Theorem  1 and  its  corollary  provides  the  bases  for  the 
proposed  algorithm  for  solving  Problem  2. 

3.  COMPUTING  OPTIMAL  SOLUTIONS  FOR  PROBLEMS  2 AND  3 


We  begin  this  section  by  developing  a method  for  determining  the  optimal 
solution  to  Problem  3.  Observe  that  the  optimal  solution  must  satisfy  the 
following  two  conditions: 


and 


dB. 

i + e c 

dN.  2 l 


0 


i=l 


c . N . 

l l 


C. 


The  second  condition  must  hold  since  each  B^(N\)  is  a strictly  decreasing 
function  of 


B. (N. ) = a.e  ^i^i , 
li  l 


Since 


•rrs«r 


- . 


where  g.  = 6d.  -a  and  f.  = 6b..  Consequently  N.  is  a linear  function 

of  C.  If  the  budget  is  incremented  by  an  amount  AC,  then  the  new  value 
of  the  stock  level  for  item  i,  N! , satisfies 


N!  = II.  + ^ . 

l l f . 

l 


The  optimal  solution  to  Problem  2 has  been  found  if  each  of  the  M. 
found  using  Equation  6 is  non-negative.  If  there  exists  an  i for  which 
fL  < 0,  then  we  may  employ  the  following  algorithn  to  find  the  optimal 
solution  to  Problem  2.  Let  I = {l,...,n}  and  represent  the  optimal 

solution  to  Problem  2. 

Step  0.  Solve  Problem  3 an  described  above  thereby  obtaining  an 
initial  value  for  h\,  iel. 

Step  1.  Set  = 0 for  all  h\  < 0 during  the  last  iteration  and 
delete  the  corresponding  i from  I.  Recompute  a and  6,  where 


c.d. 

a = l 

ini  °i 


B = 1 {c ./b. }. 

• T 11 

iel 


Step  2.  Using  Equation  6,  obtain  new  estimates  of  Ih  for  each  iel. 

If  N.  > 0 for  all  iel,  then  the  optimal  solution  has  been  found,  and 

= N.  for  all  iel  and  n}  = 0 for  all  i = l,...,n  for  which  it I.  If 
ii  l 

there  exists  some  i for  which  < 0,  return  to  Step  1. 


mm 


'mm 


It  is  clear  that  our  solution  satisfies  all  the  optimality  conditions 
for  Problem  2 except  possibly  condition  (a)  for  i£l.  However,  at  an  earlier 
iteration  (when  i was  deleted  from  I)  we  had 


1 1 N . <0 


+ Vi = °» 


where  0_  is  the  earlier  va]ue  of  Since  -r—  is  clearly  incr  easing 

^ 2 QlJ. 

x 

in  Ih,  and  6^  increases  at  each  iteration  (Theorem  1 and  its  corollary), 
condition  (a)  must  hold.  Convergence  is  guaranteed  since  n is  finite. 


4.  A COMPARISON  OF  ALTERNATIVE  SOLUTION  PROCEDURES  FOR  SOLVING  PROBLEM  .1 


In  this  section  we  briefly  review  three  algorithms  for  solving  Problem 
1 and  compare  them  to  two  algorithms  designed  to  obtain  a solution  for  Problem 
1 based  on  the  solution  to  the  aprroximating  problem , Problem  2. 

The  first  algorithm  we  will  discuss  is  the  procedure  originally  proposed 
by  Sherbrooke  [3].  It  is  a marginal  analysis  algorithm  consisting  of  two 
phases.  In  the  first  phase,  each  item  is  examined  independently.  The 
optimization  problem  solved  for  item  i in  the  first  phase  has  the  form 


minimize  T T (x  - s.  . )p(x| .T. . (s ..)) 

j=l  x>s..  ^ 13  ^ 10 

il 


subject  to 


Jo  ■«  = "*• 


s • • ~ 0,1,..., 

11 


(7) 


gfiiMMM1 


* 

i 


i 


where  is  the  total  system  stock  available  for  distribution  among  the 

depot  and  bases.  Let  Z^(N^)  represent  the  optimal  objective  function  value 
given  units  are  available  for  distribution.  Problem  7 is  solved  by 

obtaining  the  solution  to  the  h\+l  problems 


^iO5 


minimise  ll  <*  - s^  )p(xj  A (siQ)) 

3-1  X>Sij 


subject  to 


m 


l s..  = M. 
j=l  J 


Si0  ’ 


(8) 


s^j  - 0,1,..., 


and  s . . fixed 
lO 


for  = 0,1, Problem  8 can  be  solved  via  marginal  analysis.  Then 


Z.(N.)  = minimize  Z.(N.,s.  ),  where 
11  1 1 lO 


Si0 


s^q  — 0 9 • • • ) . 


The  second  phase  problem  is 


minimize  l Z . ( II . ) 


i=l 


i i 


subject  to 


t c.N.  < C 
i=l  1 1 " 


N.  = 0,1 

l 


W 


I- 

r. 

1 1 


i- 


r 
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Sherbrooke  [3]  suggests  a marginal  analysis  algorithm  be  used  to  find  a solu- 
tion to  this  knapsack  problem.  Clearly  other  procedures  could  be  employed  to 

obtain  an  optimal  solution.  In  any  case,  this  approach  requires  a substantial 


amount  of  storage  to  save  all  the  ( l.h  ) values.  For  moderately  sized 


problems — several  thousand  items — a storage  requirement  of  10  or  more 
words  may  be  needed  to  save  these  values.  Furthermore,  the  computation  time 
required  to  obtain  these  Z^(’L)  values  for  such  problems  is  very  large. 

Subsequently  Fox  and  Landi  proposed  a Lagrangian  algorithm  for  solving 
Problem  1 [2].  In  particular,  they  formulated  the  relaxed  version  of 
Problem  1 as 


m n 


min  l l l (x  - s )p(x,'X..T  (s.n)  + 6 £ l ci 


j=l  i=l  x>s 


il 


ij  ij  iO 


.s. . 


il 


j=0  i=l 


il 


(9) 


s. . = 0,1,. .. . , 


where  0 is  the  Lagrangian  multiplier.  Since  pi^blem  9 is  separable  by 
item,  its  optimal  solution  can  be  found  bv  solving  the  n individual  item 
problems 


m 


minimize  / Y (x  - s . . )p(x| 1 . .T. . (s . „ ) +0  Y c.s.. 

j = l *>s. . 13  13  ^ 10  j£0  1 13 

J il 


subject  to 


s. . = 0,1,  ... 
il 


J 


This  problem,  like  Problem  8 in  Sherbooke's  two-phase  method,  is  solved 
using  a partitioning  procedure;  that  is,  it  is  reformulated  as 


f 


(m 

0C. S.  + 1 

. 1 i=1 


;nm  • { > (x  - s. . )p(x  X.  .T. . (s  )) 

„ „ , L r 1 r it  lO 

s. . = 0,1, x>s. . J J 

-i  i] 


+ 6c.s. . : s.„  fixed} 

l i]  i0 


ixed}"'|  , 


or  equivalently  as 


minimize  Z(s^g;6) 


3i0  " 0,1,‘ 


where 


Z(Si0  6)  = 0CiSiO 


+ l min  { l (x-s..)p(:c|X  T..(si0))+  Oc.s..  : s..  = 0 ,1, . . .:si0  fi 

j=l  ij  x>s..  j j j j j 


fixed} . 


To  determine  Z(s^Q; 0),  solve  the  m base  problems 


minimize  Y (x-s..)p(:th..T..(s.„))  + ec.s... 

**  i]  1 13  ij  iu  1 ij 


s..  x>s.. 

11  11 


The  optimal  s^  is  the  smallest  non-negative  integer  for  which 


I P(x|XijTij(si0))  < 0c.. 


Problem  10  (or  Problem  11)  is  solved  for  each  item  for  a given  value  of 
6.  This  yields  a total  investment  cost  corresponding  to  0.  In  the  Fox-Landi 
approach,  the  ’optimal  value  of  0 is  selected  from  a grid  of  M equally 
spaced  values 


"o'6!*  -•  > V °- 
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The  optimal  value  of  6 is  the  0 , K e whose  corresponding  total 

K 

investment  cost  is  closest  to  C.  ^ 

Fox  and  Landi  suggest  that  their  method  is  a single  pass  method,  that  is, 
only  one  pass  through  the  item  data  base  is  necessary  to  obtain  the  optimal 
solution.  The  storage  requirement  to  effect  this  one  pass  approach  is  poten- 
tially enormous.  For  a moderately  sized  problem  having  3000  items,  20  bases 
and  M = 63,  almost  4 million  item  stock  levels  must  be  saved  plus  possibly 
millions  of  additional  item  data  elements  reflecting  fill  rates,  probability 
of  no  stockout  at  an  arbitrary  point  in  time,  expected  base  backorders,  etc. 
Furthermore,  there  may  be  no  simple  method  for  estimating  suitable  bounds  on 

the  values  of  the  multipliers  thereby  requiring  much  larger  values  of  M to 
insure  adequate  approximation  of  the  budget. 

In  the  author's  experience,  the  ability  of  Air  Force  personnel  to  estimate 
a reasonable  range  for  0 for  large  problems  is  not  good.  It  is  not  surprising 

that  it  is  difficult  for  someone  to  estimate  the  optimal  value  of  the  multiplier. 

The  data  used  in  the  model  frequently  change  in  real  situations  thereby 
causing  the  optimal  value  of  the  multiplier  to  change.  Furthermore,  changing 
the  multipliers ' magnitude  by  10  6 or  less  often  causes  the  corresponding 
total  cost  to  change  by  many  millions  of  dollars.  Consequently,  values 
of  6 have  been  used  in  some  Air  Force  applications  to  make  the  system  'fool 
proof.'  In  these  cases  60  million  or  more  item  stock  levels  would  be  needed 
to  be  stored  explicitly — plus  a considerable  amount  of  other  item  and  base 
data — to  make  the  Fox- Landi  algorithm  a truly  one  pass  method. 

On  the  other  hand,  if  their  method  is  altered  so  that  the  item  data  are 
examined  a second  time,  it  is  possible  to  eliminate  virtually  all  the  requirements 

for  secondary  storage.  In  the  first  pass,  only  the  running  total  cost  correspon- 

ding to  each  0 , Ke{0,. . . ,tl} , is  saved.  At  the  end  of  this  phase  the  ’optimal1 
K 
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multiplier  value,  0 , is  established.  The  second  phase  of  the  algorithm 
requires  a second  pass  through  the  data  base.  In  the  second  pass,  the  optimal 
stock  levels  for  each  location  are  found  for  all  items  by  resolving  Problem 
10  with  0=0. 

In  some  applications  the  Fox-Landi  one-pass  method  is  clearly  infeasible, 
that  is,  there  may  not  be  enough  peripheral  storage  capacity  to  save  all  the 
data.  If  storage  capacity  is  available,  there  is  a tradeoff  between  the  time 
and  cost  required  to  store  and  access  the  data  in  the  secondary  memory  using 
the  one-pass  method  and  the  time  and  cost  to  recompute  the  stock  levels 
using  the  second  method.  For  realistic  Air  Force  problems,  the  two-pass  method 
appears  to  be  the  only  feasible  approach  given  current  hardware  constraints  if 
M is  large  enough  to  guarantee  that  a solution  can  be  found  that  closely 
approximates  the  target  budget.. 

A third  way  to  solve  Problem  1 is  a slight  modification  of  the  Fox-Landi 

algorithm.  This  third  method,  called  the  bisection  method,  employs  a 

bisection  search  to  find  the  optimal  value  for  0.  This  procedure  requires 

initial  upper  and  lower  bounds  on  the  optimal  value  of  0-  Call  these  0^ 

and  0 , respectively.  The  bisection  method  is  as  follows: 

L» 

1.  Set  ? = (0U  + 0l)/2. 

2.  Solve  Problem  10  with  0=0  for  each  item. 

3.  If  the  total  cost  of  the  solution  obtained  in  Step  2 exceeds  C, 
then  replace  0^  with  0-  otherwise,  replace  0^  with  0. 

4.  If  a stopping  criteria  has  not  been  met  (such  as  a fixed  number  of 
iterations  or  an  error  tolerance),  return  to  Step  1 otherwise,  stop. 

The  major  drawback  to  the  bisection  approach  is  that  a separate  pass 
through  the  item  data  base  is  required  at  each  iteration  of  the  algorithm. 
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This  algorithm  performs  very  well  in  terms  of  convergence  and  in  our  experience 
virtually  always  produces  solutions  that  are  within  1/2  percent  of  the  target 
budget  using  10  bisections. 

The  closeness  of  the  solutions  to  the  target  budget  generated  by  either 
the  Fox-Landi  method  of  the  bisection  algorithm  depends  on  how  broad  a range 
of  multiplier  values  must  be  searched  for  a fixed  value  of  M or  a fixed 
number  of  bisections.  It  should  be  pointed  out  that  both  of  these  methods  only 
yield  an  approximation  to  the  optimal  multiplier  value  (assuming  one  exists). 

Of  the  methods  discussed  thus  far,  it  has  been  the  experience  of  both 
the  author  and  Fox  and  Lanai  [2]  that  the  latter  two  algorithms  dominate 
Sherbrooke’s  algorithm  in  run  times  by  an  order  of  magnitude  or  more  on  real 
problems  given  reasonable  estimates  of  upper  and  lower  bounds  for  the  Lagrangian 
multiplier.  Thus  in  the  comparisons  we  will  report,  only  these  two  Lagrangian 
methods  will  be  discussed. 

Earlier  we  described  an  approximation  method  for  estimating  the  optimal 
values  of  0 and  each  IK.  Several  options  are  open  for  implementing  this 
approximation  method.  One  way  to  implement  it  is  to  use  a two-phase  approach. 
Call  this  approach  the  First  Approximation  Method.  The  values  of  a^  and  b^ 
are  computed  in  the  first  phase  of  this  method  during  which  the  optimal  value 
of  G is  also  estimated  using  Equation  5.  In  the  second  phase,  we  solve 
Problem  10  for  each  item  using  the  estimate  of  the  optimal  0.  This  approach 
has  two  major  advantages  over  the  Fox-Landi  method: 

(a)  The  estimate  of  the  optimal  multiplier  can  be  obtained  without 
prespecifying  a range  of  values,  and  computation  time  to  obtain 
the  estimate  does  not  depend  on  the  uncertainty  of  the  multiplier 
value. 

(b)  The  computation  time  to  find  an  estimate  of  the  optimal  multiplier 
is  much  smaller. 


j 
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If  the  two-pass  version  of  the  Fox-Landi  algorithm  is  used,  the  second  phase 
of  that  method  and  the  second  phase  of  the  approximation  method  are  identical. 
The  ore-pass  version  of  the  Fox-Landi  algorithm  requires  considerably  more 
storage,  and  also  requires  more  computer  time  to  determine  the  optimal  stock 
levels  than  this  approximation  method  requires. 

The  First  Approximation  Method  also  has  the  following  advantages  over  the 
bisection  method: 

(a)  Cniy  two  passes  through  the  data  base  are  required  as  opposed  to 
seven  or  more  required  for  the  bisection  method  in  practice. 

(b)  No  stock  levels  need  to  be  saved-  in  the  bisection  method  it  is 
necessary  to  save  all  stock  levels  and  other  data  for  three 
multiplier  values. 

Another  algorithm  can  be  employed  that  directly  uses  the  results  of 
the  approximation  problem,  that  is.  Problem  2.  Call  this  approach  the  Second 
Approximation  Method.  This  algorithm  is  of  interest  in  situations  in 
which  we  only  want  to  compute  total  system  stock  for  each  item  and  are  not 
particularly  interested  in  computing  the  optimal  distribution  of  the  assets. 
Determining  the  optimal  allocation  of  a budget  among  items  is  of  primary 
importance  when  purchasing  inventory  or  making  budgetary  projections  for  spares 
for  different  systems.  In  these  cases,  distribution  decisions  are  usually  not 
that  critical. 

This  Second  Approximation  algorithm  also  consists  of  two  phases,  in  the 
first  phase  we  estimate  the  values  of  the  a^  and  b^  parameters , and  in 
the  second  phase  we  determine  total  system  stock  for  each  item  using  the 
algorithm  described  in  Section  3 and  rounding  hh  to  the  nearest  integer.  The 
algorithm  requires  one  pass  through  the  item  data  base  and  one  pass  through  an 
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item  file  consisting  of  a^,b^,  and  c^.  The  major  advantage  of  this  approach 
is  that  it  eliminates  the  stock  allocation  phase  of  both  the  Fox-Landi 
algorithm  and  the  First  Approximation  Method. 

5 . A COMPUTATIONAL  COMPARISON  OF  VARIOUS  ALGORITHMS 

The  Fox-Landi  algorithm,  bisection  algorithm,  and  the  two  approximation 

methods  have  been  coded  and  tested  on  several  sample  sets  of  data  for  the 

Air  Force's  new  F-15  fighter.  Since  all  of  the  test  yielded  the  same  general 

results,  we  will  discuss  only  two  of  them.  The  first  test  consisted  of  a 75 

item  sample  and  had  3 operating  bases.  The  flying  programs  were  very  different 

at  each  base.  In  the  second  test,  125  items  were  included  in  the  sample  with 

demands  occurring  at  5 bases.  In  the  second  test,  only  the  Fox-Landi  and  the 

two  approximation  methods  were  compared.  In  all  Fox-Landi  calculations,  a 

maximum  of  128  multiplier  values  were  examined:  ten  bisections  were  used  in 

all  applications  of  the  bisection  method.  The  run  times  stated  for  both 

approximation  algorithms  include  the  time  required  to  estimate  the  values  of 

a.  and  b..  Furthermore,  in  both  test  cases  all  stcck  levels  for  all 
i l 

relevant  multiplier  values  were  stored  in  main  memory.  Thus  the  reported  com- 
putation times,  which  include  compile  times  which  are  roughly  equal  for 
all  the  algorithms,  are  biased  in  favor  of  the  Fox-Landi  method  since  for 
larger  problems  this  type  of  storage  would  be  impossible.  Additionally,  the 
range  of  multiplier  values  considered  in  the  test  of  the  Fox-Landi  and 
bisection  methods  was  selected  after  estimating  the  optimal  multiplier  value 
usintf  the  First  Approximation  Method.  Thus  the  test  results  are  biased  in 
favor  of  them,  since  the  range  of  multiplier  values  was  much  smaller  than 
would  normally  be  the  case. 
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The  data  displayed  in  Tables  I and  II  indicate  how  well  each  approach 
approximates  a given  target  budget  for  the  two  test  data  sets.  Without  a 
doubt  the  bisection  method  produced  solutions  that  Lest  marched  the  tar  *et 
budgets  followed  in  order  by  the  Second  Approximation  Method,  the  Fox- 
Landi  method,  and  the  First  Approximation  Method.  As  mentioned  before,  the 
results  are  biased  in  favor  of  both  the  Fox-Lnndi  and  bisection  methods 
due  to  the  initialisation  of  the  range  of  multiplier  values.  From  a 
practical  viewpoint,  all  approaches  worked  acceptably  well  in  meeting  the 
target  budgets.  Furthermore,  the  stock  levels  generated  by  the  various 
approaches  were  virtually  the  same  for  similar  budgets.  Consequently,  total 
system  expected  backorders,  for  all  practical  purposes,  are  indistinguishable: 
that  is,  the  backorder  versus  investment  curves  virtually  coincide  among  these 
various  approaches.  Exact  comparison  of  computed  stock  levels  and  expected 
backorders  cannot  be  made  among  the  competing  methods  since  the  allocation 
of  the  available  budget  in  each  case  depends  on  the  way  each  algorithm 
estimates  the  Lagrangian  multiplier. 

The  area  in  which  the  methods  clearly  differ  is  in  computation  time.  The 
approximation  methods  require  substantially  less  time  than  either  the  Fox- 
Landi  method  or  the  time  consuming  bisection  method.  Other  experimentation 
has  shown  that  the  percentage  difference  in  computation  times  tends  to  be 
even  more  substantial  as  the  number  of  items  considered  increases. 

Thus  the  approximation  methods  produce  answers  that  are  as  good  as  those 
produced  by  either  the  Fox-Landi  method  or  the  bisection  method,  but  with 
less  computational  effort.  The  bisection  method  did  match  target  budgets 
slightly  better  than  the  approximation  methods.  However,  the  approximation 
algorithms  are  virtually  fool-proof.  This  is  perhaps  the  greatest  advantage 
of  the  approximation  algorithms.  The  user  does  not  have  to  specify  the  range 

__ _ __ _ J 
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of  multiplier  values  or  the  number  of  bisections  in  advance.  This  eliminates 
one  of  the  main  difficulties  associated  with  implementing  either  the  Fox-Landi 
or  bisection  algorithms.  In  view  of  these  observations,  the  approximation 
procedures  developed  here  appear  to  be  superior  for  use  on  real  problems. 


Table  2 


12 5- ITEM,  5- EASE  TEST  CASE 


Target  > 

Budget 

Total  Cost  (millions  ' 

of  $) 

Fox-Landi 

» 

Approx.  I | 

Approx . I I 

26.4 

26.7 

I 

24.8 

26.6 

27.6 

27.6 

26.2 

27.9 

28.7 

28.7 

27.6 

28.9 

29.8  i 

30.0  i 

25.5 

! 

29.8 

31.0 

31.2 

30.7 

30.8 

32.1 

i 

32.1 

32.0 

32.2 

33.2 

CO 

CO 

CO 

33.1 

, 33.1 

34.4 

34.4  j 

34.3 

34.2 

35.4 

35.5 

35.9  ; 

35.7 

36.6 

36.8 

j 

| 

37.0 

36.7 

37.8 

38.0 

38.1 

37.7 

38.9 

l o 

CO 

CO 

39.3 

39.2 

40.0 

39.9 

40.6 

1 

40.0 

41.2 

41.1 

42.1 

I 

41.3 

42.3 

42.5 

43.9 

42.4 

43.4 

43.3 

44.7 

43.7 

44.6 

44.5 

45.6 

44.2 

45.7 

46.3 

45.1 

45.9 

46.8 

hi. 2 

47.3 

46.7 

Execution 

Times 

(Seconds) 

36.98 

16.28 

4.74 

NOTE:  All  programs  were  run  on  an  IBM  370/168  using  the 

VJATFIV  compiler. 
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inventory  systems  J^jT*  Subsequently  Fox  and  Landi  ££]  proposed  a Lagrangian 
approach  for  obtaining  item  stock  levels  for  each  location.  In  this  paper  we 
develop  a method  for  estimating  the  value  of  the  optimal  Lagrangian  multiplier 
used  in  the  Fox-Landi  algorithm  present  alternative  ways  for  determining  system 
stock  levels,  and  compare  these  proposed  approaches  with  the  Fox-Landi  
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